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We compare lattice data for the short-distance part of the static energy in 2 + 1 flavor quan- 
tum chromodynamics (QCD) with perturbative calculations, up to next-to-next-to-next-to leading- 
logarithmic accuracy. We show that perturbation theory describes very well the lattice data at short 
distances, and exploit this fact to obtain a determination of the product of the lattice scale ro with 
the QCD scale Aj^jg-. With the input of the value of ro, this provides a determination of the strong 
coupling a 3 at the typical distance scale of the lattice data. We obtain a 3 (1.5GeV) = 0.326 ±0.019, 
which provides a novel determination of a s with three-loop accuracy (including resummation of the 
leading ultrasoft logarithms), and constitutes one of the few low-energy determination of a s avail- 
able. When this value is evolved to the Z-mass scale Mz, it corresponds to a s (Mz) = 0.1156tQ Q022- 

PACS numbers: 12. 38. Aw, 12.38.Bx, 12.38.Gc 



The static energy in quantum chromodynamics 
(QCD), i.e. the energy between a static quark and a 
static antiquark separated by a distance r, is a basic ob- 
ject to understand the behavior of the theory [l| and 
constitutes a fundamental ingredient in the description 
of many physical processes [2j. The short-distance part 
of the static energy can be computed using perturbative 
techniques, and it is nowadays known at next-to-next- 
to-next-to leading-logarithmic (N 3 LL) accuracy, i.e. in- 
cluding terms up to order a 4+n In" a s with n > 
(lno: s terms appear due to virtual emissions of ultrasoft 
gluons, which can change the color state of the quark- 
antiquark pair [l(J [HI)- It can also be computed on 
the lattice, and the comparison of the two approaches 
tests our ability to describe the short-distance regime of 
QCD, besides providing information on the region of va- 
lidity of the perturbative weak-coupling approach [l2[ • A 
comparison of the static energy at N 3 LL accuracy with 
quenched lattice data [l]| was presented in Ref. Q . Here 
we present lattice data for the short-distance part of the 
static energy in 2 + 1 flavor QCD and compare it with the 
perturbative calculation up to N 3 LL accuracy. This al- 
lows us to determine the strong coupling a s at three-loop 
accuracy (including resummation of the leading ultrasoft 
logarithms), in a way which is largely independent from 
the other determinations that currently enter in the world 
average [lij ]. The natural scale where our determination 
is performed corresponds to the inverse of the typical dis- 
tance where we have lattice data, i.e. around 1.5 GeV. 
Therefore, our analysis provides a determination of a s 
at a scale smaller than those entering the current world 
average 14|, and constitutes in this way an important 
ingredient to further test asymptotic freedom in QCD. 

The static energy has been calculated on the lattice 
in 2 + 1 flavor QCD using a combination of tree-level 



improved gauge action and highly-improved staggered 
quark (HISQ) action [l5| in Ref. [lq |. The strange-quark 
mass m s was fixed to its physical value, while the light- 
quark masses were chosen to be mi = m s /20. These 
correspond to the pion mass of about 160 MeV in the 
continuum limit, which is very close to the physical value. 
The calculation of the static energy was performed in a 
wide range of gauge couplings 5.9 < /3 = 10/.g 2 < 7.28. 
At each value of the gauge coupling we calculate the scale 
parameters ro and ri defined in terms of the static energy 
Eo(r) as follows 
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The values of ro and r\ were given in Ref. fl6j ] for 
each p. The above range of the gauge couplings corre- 
sponds to lattice spacing 1.909/r < a -1 < 6.991/r . 
Using the most recent value ro = 0.468 ± 0.004 fm 



[16( we get 0.805 GeV < 



< 2.947 GeV. Thus 



we can study the static energy down to distances 
r = 0.14ro or r ~ 0.065 fm. For the comparison 
with perturbation theory the most relevant data set is 
the one corresponding to lattice gauge coupling j3 = 
6.664, 6.740, 6.800, 6.880, 6.950, 7.030, 7.150, 7.280, which 
is what we will use here. The static energy can be calcu- 
lated in units of ro or r\. Since the static energy has an 
additive ultraviolet renormalization we need to normal- 
ize the results calculated at different lattice spacings to 
a common value at a certain distance. We fix the static 
energy in units of r to 0.954 at r = r [l6j]. At distances 
comparable to the lattice spacing the static energy suf- 
fers from lattice artifacts. To correct for these artifacts 
we use tree level improvement. From the lattice Coulomb 
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potential 
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d 3 k 
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we can define the improved distance 77 = (47rCx(r)) _1 
for each separation r. Here Z?oo is the tree level gluon 
propagator for the a 2 improved gauge action. The tree 
level improvement amounts to replacing r by r / 13[ . Al- 
ternatively following Ref. 
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we fit the lattice data 
at short distances to the form const~a/r + ar + a'(l/r — 
1/rj) and subtract the last term from the lattice data. 
We have found that both methods of correcting for lat- 
tice artifacts lead to the same results within errors of the 
calculations. Furthermore, the static energies calculated 
for different lattice spacings agree well with each other 
after the removal of lattice artifacts. The corrected lat- 
tice data obtained for several lattice spacings are shown 
in Fig.Q]as the points. All the lattice data seem to lie on 
a single curve even at short distances, indicating that the 
above procedure of removing the lattice artifacts works. 

As mentioned before, the static energy is known at 
N 3 LL accuracy in perturbation theory. Detailed expres- 
sions for Eg were given in Ref. 0] (and references therein) 
and will not be reproduced here. For our present analysis, 
it is only important to recall that: (i) In order to obtain a 
well behaved perturbative series, it is necessary to imple- 
ment a scheme that cancels the leading renormalon sin- 
gularity [2(| . This kind of schemes introduce dependence 
on a dimensional scale in the problem, which we denote 
as p. In particular, we implement the renormalon can- 
cellation according to the scheme described in Ref. 21 1. 
Then, the natural value for the scale p corresponds to 
the center of the range where we have lattice data (ii) At 
N 3 LL accuracy the perturbative expression depends on 
an additional constant (which was not present at lower 
levels of accuracy). This is due to the structure of the 
renormalization group equations at that order. We call 
this constant K%. It should satisfy the power counting 
condition \K%\ ~ Aj^g- (where Aj^g- is the QCD scale in 
the MS scheme), but apart from that it is unconstrained. 

We can now compare the perturbative results for the 
static energy with the lattice data. This comparison goes 
along the same lines as the quenched case in Ref. [3] , ex- 
cept that now we use rif = 3 everywhere (n f is the num- 
ber of light flavors); we also include finite strange mass 
effects at one loop, although they turn out to be negli- 
gible. We use the maximum known accuracy (four loop) 
for the running of a s everywhere (as opposed to chang- 
ing the accuracy for the running depending on the order 
we are working at), since, for the rif = 3 case, the hi- 
erarchy of scales underlying the perturbative calculation 
would not be well satisfied with the running of a s at one 
loop. The perturbative expressions depend on the value 
of the quantity roA^jg, and we will use the lattice data to 
determine it. That is, we search for the range of roAj^g- 
that is allowed by lattice data, taking into account all 
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FIG. 1. Comparison of the singlet static energy with lattice 
data (red -lighter- points) . [The comparison (and all the anal- 
ysis in the text) is done for r < 0.5ro — 0.234 fm, which is 
the region where perturbation theory is reliable. The (blue 
-darker-) points and curves for r > 0.5ro are shown just for 
illustration]. The dotted blue curve is at tree level, the dot- 
dashed magenta curve is at one loop, the long-spaced-dashed 
orange curve is at two loops, the dashed brown curve is at 
N 2 LL accuracy, the long-dashed green curve is at three loops 
plus leading ultrasoft logarithmic resummation, and the solid 



black curve is at N LL accuracy. roAj 



0.70 was used in 



all the curves. The additive constant in the perturbative ex- 
pression for the static energy is taken such that each curve 
coincides with the lattice data point at the shortest distance. 



the uncertainties involved. Then, using the value for ro 
determined in Ref. jl6} we can obtain a determination of 
a s (Mz) {Mz is the Z-boson mass). 

Under the assumption that perturbation theory by it- 
self (after canceling the leading renormalon) is enough to 
accurately describe the lattice data in the range of dis- 
tances we are considering (i.e. r < 0.5ro), a procedure to 
extract r A^g- from the comparison of the perturbative 
expressions for the static energy with lattice data was 
devised in Ref. Q . We will proceed in an analogous way 
here. The procedure exploits the fact that any value of 
p (around its natural value at the center of the range for 
which we have lattice data, i.e. p = 3.14r ( y 1 ) cancels the 
renormalon and is therefore allowed. Following Ref. 0], 
we search for a set of p values which are optimal for the 
determination of roA^jg. The procedure to do that con- 
sists of the following steps: 

1. We vary p by ±25% around its natural value p = 
3.14r f ^~ 1 , that is from p — 2.36r ~ 1 to p = 3.93TQ 1 . 

2. For each value of p and at each order in the pertur- 
bative expansion of the static energy, we perform 
a fit to the lattice data (that is we do fits at tree 
level, one loop, two loops, and three loops; in the 
last two cases with and without ultrasoft logarith- 
mic resummation). The parameter in each of these 
fits is roAj^g. 
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3. We select those p values for which the reduced x 2 
of the fits decreases when increasing the number of 
loops of the perturbative calculation. 

For the analysis, we use the fits from tree level to three 
loop plus leading ultrasoft logarithmic rcsummation ac- 
curacy. We could also use the fits at N 3 LL accuracy (i.e. 
three loops plus sub-leading ultrasoft logarithmic resum- 
mation), which would involve the additional constant K2 
(that would also need to be fitted to the data, i.e. the fits 
involve one additional parameter at this order). We find 
that, with the present lattice data, the \ 2 as a function of 
r o^MS * s ver y flat in this case, and we cannot improve our 
extraction of roAj^g- by including the fits at N 3 LL accu- 
racy in the analysis. We interpret this as the unquenched 
lattice data not being accurate enough to be sensitive to 
subleading ultrasoft logarithms (unlike the quenched one 
used in |3j), a fact that leaves room for future improve- 
ments. Therefore, we take the numbers at three loop 
plus leading ultrasoft logarithmic resummation accuracy 
as our best result, and consider, at this order, the set 
of fitted values of roA^g- for the p range obtained after 
step 3 above (we denote these values by xi). We assign 
a weight to each Xi, given by the inverse of the reduced 
X 2 of the fit. We take the weighted average of the Xi as 
our central value for the determination of roA^jg. To es- 
timate the error that we should associate to this number, 
we consider the weighted standard deviation of this set 
of values, and the difference with the weighted average 
computed using the result at the previous perturbative 
order (with the corresponding p range that one obtains 
at that order; for illustration, we show the results for 
r o^MS obtained at different levels of accuracy in Tab.H|. 
We obtain roA^ = 0.7024±0.0011±0.0665 = 0.70±0.07, 
where the first error is due to the weighted standard de- 
viation, the second to the difference with the two-loop 
result, and we summed the two errors linearly on the 
right-hand side of the equation. It is important to point 
out that the error assigned to the result must account for 
the uncertainties associated to the neglected higher-order 
terms in the perturbative expansion of the static energy; 
in that sense, assigning the difference with the result at 
the previous order as an error (as we do) is a quite con- 
servative estimate. Note that, starting at the two-loop 
level, one can decide whether to perform resummation of 
the ultrasoft logarithms or not; when assigning the error, 
we take whichever difference is larger. We also mention 
that there is an error associated to each of the Xi coming 
from the fit to the lattice data, but the error that this in- 
duces in the average can be neglected. To further assess 
the systematic errors stemming from our procedure, we 
have redone the analysis using p- value weights, obtaining 
roA^jg = 0.7022±0.0011±0.0628 = 0.70±0.06, and using 
constant weights, obtaining roAj^g = 0.7022 ± 0.0011 ± 
0.0666 = 0.70 ±0.07, which are both compatible with the 
previous analysis. In our final result we quote an error 



Accuracy 




>iAms 


tree level 


0.395 


0.397 


1 loop 


0.848 


0.862 


2 loop 


0.636 


0.654 


N 2 LL 


0.756 


0.783 


3 loop 


0.690 


0.701 


3 loop + us. res. 


0.702 


0.715 



TABLE I. Values of roA^jg- obtained at different levels of accu- 
racy. The second column shows the results obtained using the 
static energy normalized in units of the scale ro. The third col- 
umn shows the results obtained with the static energy normal- 
ized in units of the scale ri and then (for easier comparison) 
transformed to ro units, using the factor ro/ri = 1.508±0.005 
[l6| . "N 2 LL" stands for next-to-next-to leading-logarithmic 
(i.e. two loop plus leading ultrasoft logarithmic resummation) 
and "3 loop + us. res." stands for three loop plus leading 
ultrasoft logarithmic resummation. 



that covers the whole range spanned by the three analy- 
ses. As an additional cross-check of the result, we have 
redone the analysis with the static energy normalized in 
units of the scale r\ — 0.3106 ± 0.0020 fm (rather than 
ro); these numbers are presented in the third column of 
Tab. [H and are consistent with our previous results. Fi- 
nally, we point out that in some cases the \ 2 as a function 
of roA-jgg at next-to-next-to leading- logarithmic (N 2 LL) 
accuracy develops a second local minimum for larger val- 
ues of ro Aj^g. To discern which minimum should be taken 
as the physical result when this happens, we have redone 
the fits using smaller r ranges. We found that the posi- 
tion of the second minimum is not stable, while the po- 
sition of the first one is. Furthermore the minima of the 
X 2 from lower orders in perturbation theory come closer 
to the first of the two minima at N 2 LL, when decreasing 
the r range we use. In view of the above, when a second 
minimum develops, we keep the first one, which is stable 
and preferred by lower perturbative orders. 

According to the discussion in the previous paragraph, 
our final result reads 

r A m = 0.70 ± 0.07, (3) 

which using ro = 0.468 ± 0.004 fm [l6| gives 

a s (p = 1.5GeV,n/ = 3) = 0.326 ± 0.019, (4) 

the uncertainty in ro is negligible in the final error above. 
When we evolve Eq. ^ to the scale Mz we obtain 



a s (M z ,n f = 5) = 0.1156 



+0.0021 

0.0022' 



(5) 



where we have used the Mathematica package RunDec 
[22I ] to obtain the above number (4 loop running, with 
the charm quark mass equal to 1.6 GeV and the bottom 
quark mass equal to 4.7 GeV). We mention that the final 
result employing the static energy normalized in units of 
n is a s (M z ) = 0.1160Iq;oo22; which is compatible with 
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our result in Eq. §5§ and further shows its robustness. 
Figure [T] shows a comparison of the perturbative expres- 
sions for the static energy with lattice data using our 
result ToAjvfs = m Eq. © (p is set at the natural 
value, p — 3.14r _1 ). We can see that the perturbative 
series converges, approaches the lattice data, and repro- 
duces it very well at N 3 LL accuracy (the constant K2 
that appears at N 3 LL is fixed by a fit to the lattice data, 
which gives tqK2 — —2.3, fulfilling the power counting 
condition) . Note also that our results for a s are not sen- 
sitive to the specific value of r that we consider as the 
upper limit where perturbation theory is reliable, as it is 
manifest from Fig. [TJ 

Our result in Eq. ([5]) constitutes a novel determina- 
tion of a s (since it is largely independent of the other 
available determinations), that stems from a perturba- 
tive calculation of the QCD static energy at three loop 
plus leading ultrasoft logarithmic resummation accuracy. 
With respect to the other determinations currently en- 
tering the world average it represents the one at lowest 
energy. The lowest-energy determination so far was that 
from the r system, performed at m T — 1.78 GeV. Our 
result is therefore an important new ingredient to test 
the running of a s . 

Other recent determinations of a s , that also employ 
comparisons with lattice data, include Refs. [23], [24| 
where several observables related to Wilson loops (but 
not the static energy) are used, Refs. 2J, |25j which em- 
ploy moments of heavy quark correlators, Ref. |26l| that 
uses the vacuum polarization function, Ref. [271 ] which 
uses the so-called Schrodinger functional scheme (albeit 
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employing rather high pion masses), and Ref. [28[ that 
employs the ghost-gluon coupling; they deliver numbers 
that are mostly compatible with our result, although 
with central values a bit higher than ours. We also men- 
tion that comparisons of perturbative calculations for the 
static energy with lattice data in QCD with n f = 2 fla- 

II S3. 



vors have been presented recently in Refs. 

Let us also point out that the comparison of the pertur- 
bative result with lattice data, shown in Fig.[TJ is interest- 
ing in itself, since the static energy constitutes a basic in- 
gredient in the description of many physical processes 0] . 
As an example, it is relevant for the study of quarkonium 
production in heavy-ion collisions: A direct lattice cal- 
culation of the quarkonium spectral functions is known 
to be difficult 3l|. The only viable option could be cal- 
culating quarkonium spectral functions within an effec- 
tive field theory framework, as potential Non-Relativistic 
QCD (pNRQCD) provides [H[33|]. To access the valid- 
ity of a weak coupling pNRQCD approach at non-zero 
temperature, one eventually will have to compare weak- 
coupling calculations of the static quark- antiquark cor- 
relators with the corresponding lattice calculations. The 
comparison of the static energy at zero temperature to 
the perturbative results, in 2+1 flavor QCD with physi- 
cal quark masses (as we have provided), is an important 
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FIG. 2. Determinations of a s that enter in the world average 
[l4| (blue -darker- points) compared with our result (most- 
left red -lighter- point), as a function of the energy scale Q. 
The band is the world-average value of a s (Mz), evolved with 
4-loop accuracy. 



first step in this direction. 

In summary, we have shown that perturbation theory 
(after canceling the leading renormalon singularity) can 
describe the short-distance part of the QCD static energy 
computed on the lattice, see Fig. Q] this is the first time 
that this is done for QCD with rij = 2+1 dynamical 
quarks. Exploiting this fact, we have obtained the range 
of roAj^g that is allowed by lattice data. Using the value 
of ro as an additional input, this provides a determination 
of a s . We obtained a s (1.5GeV) = 0.326 ± 0.019, which 
represents one of the few low-energy determinations of a s 
available and, when evolved to the scale Mz, corresponds 



to a s (M z ) = 0.1156 



+0.0021 

-0.0022- 



A comparison of our result 



with the determinations of ot s that currently enter in the 



world average 14J is shown in Fig. [21 
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